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Abstract: 

Improving the accuracy of the galactic cosmic ray (GCR) environment and transport 
models is an important goal in preparing for studies of the projected risks and the 
efficiency of potential mitigations methods for space exploration. In this paper we 
consider the effects of the isotopic composition of the primary cosmic rays and the 
isotopic dependence of nuclear fragmentation cross sections on GCR transport models. 
Measurements are used to describe the isotopic composition of the GCR including their 
modulation throughout the solar cycle. The quantum multiple-scattering approach to 
nuclear fragmentation (QMSFRG) is used as the data base generator in order to 
accurately describe the odd-even effect in fragment production. Using the Badhwar and 
O’Neill GCR model, the QMSFRG model and the HZETRN transport code, the effects of 
the isotopic dependence of the primary GCR composition and on fragment production for 
transport problems is described for a complete GCR isotopic-grid. The principle finding 
of this study is that large errors (+ 100%) will occur in the mass-flux spectra when 
comparing the complete isotopic-grid (141 ions) to a reduced isotopic-grid (59 ions), 
however less significant errors (<30%) occur in the elemental-flux spectra. Because the 
full isotopic-grid is readily handled on small computer work-stations, it is recommended 
that they be used for future GCR studies. 


INTRODUCTION 


In the description of the transport of the galactic cosmic rays (GCR) in shielding 
materials or tissue, a common approximation is to consider only the elemental 
composition of the primary GCR and a reduced isotopic-grid for the secondary nuclei 
produced in nuclear fragmentation interactions. In this paper, we analyze the role of the 
isotopic dependence of the GCR primary composition and nuclear fragmentation in 
predicting the flux of the GCR behind arbitrary shielding configurations. NASA has set a 
strategic goal of determining the GCR flux spectra to within +25% accuracy (Annon. 
1998). Our study is an important milestone in achieving this goal, since for the first-time 
a full isotopic-grid has been achieved in a GCR transport model and we document the 
error inherent in former approaches. Also, for applications that will consider radioactive 
isotopes produced in the atmosphere or shielding, our study provides a useful tool to 
perform such analyses. 

Historically the HZETRN code grew from a 29-isotopic grid used in the 1980’s and early 
1990’s (Wilson, et al., 1991) to an extension to 32 isotopic-grid in order to include all 
light ions (Cucinotta, 1993). Because of the limitations of random access memory present 
in the computer workstations of the early 1 990’s, sensitivity studies were made for mono- 
energetic ion beams to study the minimum number of isotopes for convergence resulting 
in the use of 59-isotopic grid (Kim, et al., 1994), and all GCR studies since 1994 have 
used the 59-isotopic grid (Shinn, et al., 1994). However, there are several reasons to re- 
consider the use of the full isotopic-grid for GCR transport problems. First, the isotopic 
dependence of the primary GCR has not been considered in the past and may lead to 
errors in the description of both primary ion attenuation and secondary particle 
production. Secondly, the studies of Kim et al. (1994) used the NUCFRG2 model of 
fragmentation (Wilson et al., 1994), which does not provide a correct description of the 
even-odd effect observed in fragment production and of the projectile isospin dependence 
observed experimentally (Knott, et al., 1996, 1997; Zeitlin et al., 2001). Thirdly, the 
sensitivity studies made by Kim et al. 1994 used a “calibration” of the isotopic-grid to 
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j6 Fe beams, however a larger isotopic grid occurs when all GCR projectile nuclei are 
considered. Fourth, the error in the range-energy and stopping powers that results from 
the use of a reduced isotopic grid, although expected to be small for large mass number, 
A»l, is an unnecessary one for transport calculations. Finally, the improved 
computational speed and memory available on current small computer workstations, 
makes the inclusion of a full isotopic-grid in the HZETRN code to be straightforward 
problem at this time. 


GCR ISOTOPIC COMPOSITION 

Theoretical models and upper atmosphere and satellite measurements of the GCR have 
long considered the isotopic composition of the GCR in order to understand the source of 
these particles and their modification through transport through space (Parker, 1965). 
Theoretical models consider the transport of the GCR through the inter-planetary medium 
in order to extract the primary nuclear composition at stellar sources (Parker, 1965, 
Webber, et al., 1990a). Experimental studies have included measurements at the top of 
the Earth’s atmosphere and on the Pioneer I and II spacecraft. A survey of such data 
(Hesse et al., 1991; Lukasiak, et al., 1993; Webber et al., 1985, 1990; Wiedenback et al, 
1981, 1985) was made with the results shown in Table I. For this compilation, we note 
that since secondary fragment production is modulated by the time spent in the inter- 
planetary medium, the isotopic fraction is dependent on position in the solar cycle. In 
Table I we have used data on isotopic fractions near solar maximum where the isotopic 
fraction for nuclei produced within the heliosphere are expected to be at a maximum 
(Hesse et al., 1991; Lukasiak, et al., 1993; Webber, et al., 1985, 1990; Wiedenback, et al, 
1981, 1985). We then used estimates of the source abundances and the modulation 
parameter to estimate the isotopic fraction as a function of position in the solar cycle. 
Based on the results of Lukasiak et al. (1993), we use the following empirical formula to 
describe the dependence of the isotopic ratio (near-Earth) on the solar cycle as described 
using the modulation parameter <f>(MV) 


(1) f(A,Z) = / Soww (^,Z) + [exp(^)-l] 


where y is found by interpolating from the source and near solar maximum values. The 
isotopic fraction is than sealed to the Badhwar and O'Neill model spectra (1992) for the 
most abundant isotope of a given charge to obtain the primary isotopic energy spectra. 
For the Z=7 and Z=2 ions we use the empirical formula of Cucinotta et al. (1993) to 
estimate the primary (near-Earth) deuteron and helion spectra. Figure 1 shows results for 
the empirical scaling on the solar modulation parameter described by eq.(l) for l3 C and 
15 N. 

ISOTOPIC EFFECTS IN GCR TRANSPORT 

The Boltzmann transport equation within the straight-ahead and continuos slowing down 
approximations and assuming velocity conservative fragmentation events, describes the 
heavy ion flux, </>/E,x) of an ion j with mass number Aj and charge number Z, , energy E 
(in unites of MeV/u) at shielding depth x (in units of g/cm 2 ) and is written (Wilson, et al. 
1991) 

( 3 ) 

ox S j \E) oE k 

where o)(E) is the energy-dependent absorption cross section (cm' 1 ) and qj*(E) is the 
fragmentation cross section for producing an ion j from k. Because eq. (3) is a coupled 
integro-differential equation for the many GCR primary and secondary nuclei, required 
memory allocations will increase rapidly as the number of ion species is increased. 
However, such practical limitations no longer exist on small computer workstations and 
an unrestricted isotopic-grid can presently be implemented with no memory or storage 
problems. When using a reduced isotopic grid, fragments not contained in the grid are 
assigned to a nearby mass of identical charge number. This introduces two types of 
errors: First the range-energy and stopping powers are altered from their true values. This 
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error is expected to be small at high-energies for A»1 since here ion ranges are well 


fragmentation cross sections. Here because of the re-assignment of the mass number 
when using a reduced grid, an error is introduced by the change in neutron number from 
its true values. Such errors arise from the atomic transport properties, from the isotopic- 
spin dependence of fragment spectra, and in the spectra of high-energy neutrons 
produced in fragmentation events. This latter error is expected to be small because high- 
energy neutrons are largely produced by light-particle (n, p, d, t, h, and a) interactions on 
target nuclei especially for larger target atoms (A>10), but may be non-negligible for 
light target atoms (Cucinotta, et al. 1998a). Since materials with high-hydrogen content 
are known to be the optimal shielding materials, the changes in the neutron flux due to 
the used of a full isotopic-grid should be considered. 

QUANTUM ABRASION-ABLATION 

The quantum multiple scattering description of the fragmentation (QMSFRG) has been 
quite successful in describing heavy ion fragmentation data (Cucinotta, et al. 1994, 1997, 

1 998b). The prior work considered the derivation of the pre-fragment excitation spectrum 
following nucleon or alpha particle abrasion and showed that the excitation spectrum 
could be represented in terms of an impact parameter dependent convolution of the pre- 
fragment excitation response for a transition of the pre-fragment core from state n to n 
and the project fireball response, 


where b (b’) is the impact parameter, and q the momentum transfer. The abrasion 
response is defined as the interaction of the projectile fireball with the target after 
performing closure over the final fireball states 


described by A/Zf scaling factors. The second type of error occurs in the absorption and 
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(5) A 


f 

(2^r) 3 


< R'\ Q RT \b') \RxR\ Q RT (b) | R' > S(E,-E f ) 


where the Ort represent the fireball-target profile operator, and Jcr the projectile fireball 
momentum vector. The abrasion-response represents a complicated many-body operator. 
The pre-fragment excitation is described in terms of the transition matrix 

(6) P ny = <F* n \Q F . T + (b')\F*,xF\.\Q p . T (b)\F* n > 

where the matrix elements for the pre-fragment excitation are evaluated over the many- 
body profile operators. Oft ■ In the present model a convolution approach is used to derive 
the mutli-knockout spectrum using the single-fragmentation term (Cucinotta and Dubey, 
1994). 

The de-excitation of the pre-fragments in the QMSFRG model is described in a stochastic 
process using a Master equation for nuclear de-excitation by particle emission (Cucinotta 
et al., 1997). Let f (E, t) be the probability of finding the nuclei b at time t with excitation 
energy Eb and Pk (E) be the probability that the nuclei, b will emit ion k with energy E, 
then the Master equation is 

(7) df '' {E , M) = (ElDP^E) 

at / k 

In eq. (7) the first-term on the right corresponds to gains by decays a -> b+j and the 
second terms from losses dues to decays b -> c+k where the j (or k) are light-particle 
emissions ( n , p, d, t, h, or a). Equation (7) is solved by iteration up to medium excitation 
energies (<150 MeV) and by approximation for high excitation energies (Cucinotta et al. 
1997, 1998b). An important feature of the solution is the correct description of the level- 
density including the description of nuclear-shell effects at low-excitation energies, and 
the use of measured values for the nuclear masses. The fragmentation cross section is 
then evaluated from equations (4)-(7) as 
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(8) a F (A F ,Z F ) = X \ dE F> ^r~f(A F .,Z F . >A f ,Z f ) 

a f ',z,.' dE F , 

where / (A F . , Z F , » A F , Z F ) is the solution to equation (7). 

The QMSFRG theory reduces to the OPTFRG model (Townsend et al., 1986) when 
energy conservation is ignored and the closure approximations on the pre-fragment and 
fireball states are made, and to the NUCFRG2 model (Wilson, et al. 1995) when the 
optical operators are expressed as volume overlaps functions. In both the OPTFRG and 
NUCFRG2 models all information on the pre-fragment excitation spectrum is lost and 
thus must be introduced in an ad-hoc manner independent of the collision model used to 
describe mass removal. 


RESULTS AND DISCUSSION 

We first illustrate the accuracy of the QMSFRG model and the effects of isospin on 
fragmentation cross sections. Figures 2-9 show comparisons of the model to experimental 
data for the elemental distributions of fragments for several nuclei of similar mass 
number. The isospin, T : =0 nuclei display large odd-even effects, which are reduced for 
the T : * 0 nuclei. The odd-even effects are present for all target nuclei, however are 
reduced for hydrogen targets due to the small abrasion probability for large mass removal 
on hydrogen. One of the errors that can be seen in transporting ions using a reduced 
mass-grid by comparing fragmentation cross sections for nearby projectiles where large 
differences in many of the production cross sections occur for neighboring projectile 
nuclei. The model accurately reproduces the effects observed in the experiments. Figures 
1 0 and 1 1 show comparison of QMSFRG results to experiments for the isotopic 
distribution of fragments for 4(, Ar and ^Fe projectiles. The results show good agreement 
between theory and experiments and indicate the larger number of isotopes that are 
produced during nuclear fragmentation. 
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Figure 12 shows results from the HZETRN code at solar minimum ( <p=428 MV) behind 
5 g/cm 2 of aluminum shielding. Comparison of the mass-flux spectra for a 59-isotope 
grid and the full 141 -isotope grid are shown. The 141 -isotope grid was developed by 
considering the fragmentation cross sections for a large number of GCR primary nuclei 
and dominant fragments in several materials. In Figure 12 we have scaled the flux by the 
charge-squared of the ion as a measure of the ionization power of each mass group. Large 
differences are seen for many nuclei, especially odd-mass nuclei and the non-abundant 
even-mass nuclei. Figure 13 shows the percent error resulting from the use of the reduced 
isotopic-grid for shielding depths of 5 and 20 g/cm 2 of aluminum. Errors greater than 
1 00% are seen for many nuclei, however in most cases such large errors only occur for 
the less abundant nuclei. Biological effects are expected to be approximately proportional 
to Z 2 and the elemental-flux distribution may be a sufficient test of transport models for 
supporting exploration studies. In Figure 14 we show a similar comparison to that of 
Figure 13, however here for the elemental-flux distribution. The errors are indeed less 
substantial than those of the mass-flux distribution, yet are larger than 10% in many 
cases. Similar comparisons are near solar maximum conditions (<jf=1000 MV) are shown 
in Figures 15 and 16. The errors at solar maximum are slightly larger than at solar 
minimum because of the larger buildup of “secondary” nuclei in the heliosphere for these 
conditions. 

A major part of developing the description of GCR transport are theoretical models and 
experimental data on the nuclear interactions and propagation of heavy ions (Z>2). Over 
the last 30 years such descriptions have improved dramatically with the development of 
an accurate free space GCR model (Badhwar and O’Neill, 1992), the HZETRN code 
(Wilson, et al., 1990), the measurement of a significant number of fragmentation cross 
sections (for e.g., Brechtmann and Heinrech, 1988; Knott, et al. 1996, 1997; and Zeitlin 
et al., 1997, 2001, 2002), and the development of an accurate nuclear fragmentation 
model (Cucinotta et al., 1998a). Laboratory (Shimmerling, et al, 1989), and spaceflight 
(Badhwar and Cucinotta, 2000) validation data have also become available during this 
time period. The implementation of heavy transport models has progressed from models 
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that lacked unitarity (Letaw, et al. 1983) to the current fully energy-dependent models 
with accurate absorption cross sections (Shinn, et al. 1994; Wilson, et al., 1993; 
Cucinotta, 1993). Future work may still be required for light-particle transport («, p, d, t, 
h, a, and mesons and their decays), including production cross sections and studies of the 
need for 3-dimensional transport codes since angular deflections are more important than 
for heavy particles. However, the heavy ion problem is in much better shape with many 
of the remaining task ones of implementation. One exception may be further 
improvements in fragmentation cross sections and laboratory validation for the Z=J to 5 
nuclei produced from the heavier projectile nuclei (Z>10). The present paper addressed 
two implementation tasks, the use of a free-space GCR model, which includes the 
isotopic composition of the primaries, and the extension of the HZETRN code to a full 
isotopic-grid. Because the use of a reduced-grid leads to error and there are no practical 
limitations in using the full isotopic-grid at this time, we recommend that the full 
isotopic-grid be used when initiatives to design space exploration vehicles begin. Future 
tasks that remain are to implement physical models (Webber, et al. 1 990a) of the GCR 
isotopic environment and to continue to refine the QMSFRG model including 
comparisons to new fragmentation data as they become available. Such tasks are being 
considered by the present author and will be reported elsewhere. 
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Figure Captions 


Fig. 1 . Parametric model for describing the change in isotopic composition versus the 
solar modulation parameter, <p(MV). 

Fig. 2. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 20 Ne on 12 C interactions at 0.6 GeV/u. Experimental data from Zeitlin 
et al. (2001). 

Fig. 3. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 22 Ne on 'H interactions at 0.894 GeV/u. Experimental data from Knott 
et al. (1996). 

Fig. 4. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 26 Mg on 'H interactions at 0.576 GeV/u. Experimental data from Knott 
et al. (1996). 

Fig. 5. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 28 Si on l2 C interactions at 0.6 GeV/u. Experimental data from Zeitlin 
et al. (2002). 

Fig. 6. Comparisons of the QMSFRG model to experiment for the elemental distribution 
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of fragments from S on Al interactions at 1 .2 GeV/u. Experimental data from 
Brechtmann and Heinrich (1988). 

Fig. 7. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from j6 Ar on 'H interactions at 0.765 GeV/u. Experimental data from Knott 
et al. (1996). 

Fig. 8. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 40 Ar on 'H interactions at 0.352 GeV/u. Experimental data from Knott 
etal. (1996). 

Fig. 9. Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 40 Ca on 12 C interactions at 0.763 GeV/u. Experimental data from Knott 
et al. (1996). 

Fig. 10. Comparisons of the QMSFRG model to experiment for the isotopic distribution 
of fragments from 40 Ar on l2 C interactions at 0.6 GeV/u. Experimental data from Webber 
et al. (1990b). 

Fig. 1 1 . Comparisons of the QMSFRG model to experiment for the elemental distribution 
of fragments from 56 Fe on l2 C interactions at 0.6 GeV/u. Experimental data from Webber 
et al. (1990b). 


t 


Fig. 12. Comparisons of results from the HZETRN Code for the mass flux distribution 
behind 5 g/cnr of aluminum shielding for solar minimum conditions comparing transport 
with a reduced 59-isotope grid to transport with a full 141 -isotope grid. 

Fig. 13. Comparisons of the error that results from the HZETRN Code for the mass flux 
distribution near solar minimum when using a reduced 59-isotopc grid compared to 
transport with a full 141 -isotope grid. 

Fig. 14. Comparisons of the error that results from the HZETRN Code for the elemental 
flux distribution near solar minimum when using a reduced 59-isotope grid compared to 
transport with a full 141 -isotope grid. 

Fig. 15. Comparisons of the error that results from the HZETRN Code for the mass flux 
distribution near solar maximum when using a reduced 59-isotope grid compared to 
transport with a full 141 -isotope grid. 

Fig. 16. Comparisons of the error that results from the HZETRN Code for the elemental 
flux distribution near solar maximum when using a reduced 59-isotope grid compared to 
transport with a full 141-isotope grid. 
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Table la. Isotopic Composition of GCR Elements Z=5-13. 


Isotope 

Near-Earth Fraction 

Source Fraction 


Z=5 


h> b 

0.31 

0.2 

"b 

0.69 

0.8 


Z=6 


,2 C 

0.924 

0.999 

n c 

0.076 

0.001 


Z= 7 


14 n 

0.427 

0.782 

,5 n 

0.573 

0.218 


Z=8 


16 o 

0.946 

0.985 

l7 0 

0.0274 

0.008 

,8 0 

0.0265 

0.007 


Z=10 


20 Ne 

0.547 

0.681 

2l Ne 

0.106 

0.0 

22 Ne 

0.347 

0.319 


Z=12 


24 Mg 

0.648 

0.736 

25 Mg 

0.175 

0.136 

26 Mg 

0.177 

0.128 


Z=13 


26 A 1 

0.02 

0.0 

27 A1 

0.98 

1.0 
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Table lb. Isotopic Composition of GCR Elements Z=14-26. 


Isotope 

Near-Earth Fraction 

Source Fraction 


Z=14 


28 Si 

0.840 

0.902 

29 Si 

0.081 

0.054 

30 Si 

0.079 

0.044 


Z=16 


32 S 

0.698 

0.962 

33 s 

0.145 

0.019 

34 s 

0.157 

0.019 


Z=17 


35 C1 

0.523 

1.0 

36 C1 

0.414 

0.0 

37 CI 

0.261 

0.0 


Z=18 


36 Ar 

0.637 

1.0 

37 Ar 

0.314 

0.0 

38 Ar 

0.026 

0.0 

40 Ar 

0.026 

0.0 


Z=20 


40 Ca 

0.4 

1.0 

41 Ca 

0.2 

0.0 

42 Ca 

0.2 

0.0 

43 Ca 

0.2 

0.0 

44 Ca 

0.2 

0.0 


Z=26 


54 Fe 

0.076 

0.055 

55 Fe 

0.084 

0.078 

56 Fe 

0.763 

0.792 

57 Fe 

0.076 

0.075 
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Figure 1. 
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Figure 2. 


20 Ne (0.6 GeV/u)+ 12 C -> Z F 
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Figure 3. 


22 Ne (894 MeV/u)+ H -> 2 
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Figure 4. 
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Figure 5. 


28 Si (0.6 GeV/u) + 12 C -> Z F 
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Figure 6. 
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Figure 7. 


36 A r (765 MeV/u) + H -> Z F 
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Figure 8. 


40 Ar (352 MeV/u)+ H -> Z F 
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Figure 9. 


40 Ca (763 MeV/u)+ H --> Z F 



Figure 10. 





40 Ar+ 12 C (0.6 GeV/u) 
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Figure 11. 


56 Fe+ 12 C (0.6 GeV/u) 
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Figure 12. 
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Figure 13. 
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Figure 14. 
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Figure 15. 
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